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Abstract 



We derive the coordinate space wave function for the inhomogeneous six- vertex model 
^ . from the Algebraic Bethe Ansatz. The result is in agreement with the result first obtained 

long time ago by Yang and Gaudin in the context of the problem of one-dimensional 
fermions with 5- function interaction. 



Ph' 1. Introduction. 

, The solution of the six-vertex model with an arbitrary inhomogeneity parameters was an 

important step in the theory of integrable systems. In the rational case the solution was given 
independently by Yang [1] and Gaudin [2] using the generalization of the coordinate Bethe 
, ansatz (for a review see for example [3]). Later the solution was simplified drastically in the 
CN ! framework of the Algebraic Bethe Ansatz method (for example see [4]). However, since up to 
^ \ now the form of the coordinate space wave function of the eigenstates of the transfer matrix in 
^ \ the framework of the Algebraic Bethe Ansatz was not investigated, the connection of this two 



approaches remains obscure. The goal of the present letter is to fill this gap and derive the 
O ! coordinate space wave function of the Bethe eigenstate given by the Algebraic Bethe Ansatz. 
We find the complete agreement with the results of ref's [1], [2]. 

The solution of this problem is achieved with the help of the so called Factorizing operator. 
^ \ This operator introduced in ref.[5] plays an important role in the solution of Quantum Inverse 
5^ \ Scattering problem [6] , calculations of the correlation functions in spin chains [7] and scalar 
products [8]. The calculation of the wave function which looks hopeless in terms of the usual 
Algebraic Bethe Ansatz operators, is quite simple in the F-basis for the operators in the 
auxiliary space given by the tensor product of M (the number of up-spins) spin 1/2 auxiliary 
spaces. In Section 2 we introduce the notations and briefly discuss the definition and the 
properties of the Factorizing operator. We calculate the wave function of Bethe eigenstate in 
Section 3. Finally in Section 3 we present the conclusion. 

2. Algebraic Bethe Ansatz and the F-basis. 

We consider in this letter the transfer matrices corresponding to both rational and trigono- 
metric regimes of the six-vertex model. In the present section we diagonalize the operator A 
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and introduce the factorizing operator. Let us fix the notations: the normahzation of basic 
S - matrix, the definition of monodromy matrix and write down the Bethe Ansatz equations. 
For the rational case the S- matrix has the form 5*12 (ti, ^2) = — ^2 + ^-Pi2, where P12 is the 
permutation operator. In general trigonometric case it can be written as 



'S'l2(^l, ^2) 
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One can choose the normalization a{t) = 1 so that the functions b{t) and c{t) become 
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where = t for the rational case and (f)(t) = sin(t) for the trigonometric case. With this 
normalization the S'-matrix satisfies the unitarity condition Si2{ti,t2)S2i{t2,ti) = 1. The 
monodromy matrix is defined as 

To{t,m = Sio{Cl,t)S2o{C2,t)...SLo{CN,t), 

where arc the inhomogeneity parameters and L is the length of the lattice. We define the 
operator entries in the auxiliary space (0) as follows: 



/ A{t) Bit) \ _ 



a,/3 = (l,2) = (t;;). 



We denote throughout the paper (t; i) = (1; 0) so that the pseudovacuum (quantum reference 
state) |0) = 1(00. ..0}i). The triangle relation (Yang-Baxter equation) reads: 



'S'l2'S'l3'S'23 — 'S'23'S'l3'S'l2) RoO'ToTq' — Tq'TqRqqi] R 



00' 



'S'o'o- 



The action of the operators on the pseudovacuum is: A{t)\0) = a{t)\0) {a{t) = Yia'^i^a — t)): 
D{t)\0) = |0), C{t)\0) ~ 0. The Bethe Ansatz equations for the eigenstate of the Hamiltonian 
n^i-B(?i)|0) corresponding eigenvalue of the transfer - matrix Z{t) — A{t) + D{t) 

are 

M M 

(^(Qi) = n c(?a - qi){c{qi - ga))"\ Ht, {qa}) = a{t) c~^{Qa - ^) + n - Qa), 



a=l 



a=l 



where Qa are the solution of Bethe Ansatz equations. 

In the present paper we will use the monodromy matrix in the F -basis, the basis obtained 
with the help of the factorizing operator F introduced in ref.[5]. One can construct the 
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operator F = -Fi...^ which diagonahzes the operator A{t) [5], [6], [8] {A^[t) = F ^A[t)F). The 
diagonal operator A^{t) has the following form: 

A'\t) = l[m-t){l-ni)+n,). (1) 
1=1 

where we denote by rij the operator of the number of particles (hard-core bosons, corresponding 
to the up-spin) at the site i. Let us briefly mention some of the properties of the operator F. 
The explicit form of the operator F is 

F12...N = F1F2 ...Fl, Fi^{l- Hi) + Tifii, (2) 

where -hi is the operator of the number of particles (spin up) at the site i and the operator Tn 
is given by the equation 

Tn — Sn+l,nSn+2,n ■ ■ ■ ^Ln- 

One can obtain the following formulas for the matrix elements of the operator F [8] in the 
following form: 

F{m}{n} = {{m}\B{U)mn,)---B{CnMm, 

where the sets of coordinates {m} and {n} label the positions of the occupied sites. The similar 
expression can be obtained for the inverse operator F~^. Apart from diagonalizing the operator 
A{t), the operator F is the factorizing operator [5] in the following sense. For any permutation 
of indices a & Sl {Sj^ - is the group of permutations) we have the equation F = F'^R^ , where 
Fi-2..L = Ffjifj2..aL (including the permutation of the inhomogeneity parameters ^i) and Rl l is 
the operator constructed from the S- matrices defined in such a way that for the permutation 
of the monodromy matrix Tq = To^^i„2..(tL wc have = {R'^)~^ToR'^ . For the particular 
permutation cr({n}) such that crl = rii, . . . aM = Um (^1 < n2 < • • • < Um) the factorization 
condition is represented as F(F^«"»)-i = T„i-T„m- To prove the factorizing property of 
the operator (2) it is sufficient to consider only one particular permutation, for example, the 
permutation {i, i + since all the other can be obtained as a superposition of these ones for 
different i. One can show, that F — S'j+i^iF('''+^\ which evidently proves the factorization 
property. 

The matrix elements of the operators B{t) and C{t) in the F - basis: B^(t) — F~^B(t)F 
(and the same for C{t)) have the following form 

B^{t) = E^l Kii - ^) n - tWiik - i^r\l - nu) + nu) . (3) 

C^it) = ^ar b{^, - t) n (5(a - m - n,) + (c(e. - ik))-W) • (4) 

The operators (3) and (4) are quasilocal i.e. they describe flipping of the spin on a single site 
with the amplitude depending on the positions of the up-spins on the other sites of the chain. 
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The operator D^{t) can be found, for example, using the quantum determinant relation and 
has a (quasi)bilocal form. 

For our calculations performed in the next Section it will be useful to define the operators 
Bi{t), i — 1...L creating the particle (spin up) at the site i defined as B^(t) = J2iBi{t). 
According to Eq.(3) we have explicitly: 

Bi{t) = a} - t) n (cfe - - - rik) + n^) . (5) 

It is important that the operators Bi {t) have a very simple commutational relations with the 
operator A^{t') (in contrast to the well known commutational relations of B{t) and A{t')). In 
fact we have: 

B,{t)A^(t') = c(^, - t')A^{t')B,{t). (6) 

Equations (6) will be used later in the next Section. 

Let us mention here that the operators Bi — Bi{0) have the following commutational 
relations: 

From the equation (7) it easy to see that the matrices Bi {Bi = BiA^) and A^ taken in 
the auxiliary space corresponding to the lattice of M sites Oi, . . . Om with the corresponding 
spectral parameters give the explicit realization of the matrices entering the so called Matrix 
Product Ansatz for the XXZ spin chain [9] . 

3. Coordinate space wave function. 

We have to calculate the coordinate-space wave function for the inhomogeneous six-vertex 
model which defines the eigenstate |0) of the transfer- matrix according to the equation 

1^) = Yl 'iPixi, ■ ■ ■Xm)\xi, . . .Xm), 

XI, ...Xm 

where the sum is over the configurations of particles (up-spins) Xi ^ Xj for i ^ j, and 
\xi, . . . Xm) is the state with the occupied sites 1 < Xi < L. The wave function is given by the 
following scalar product: 

ip{xi, . ..Xm) = {xi, . . .XM\B{qi)B{q2) . . . 5(gM)|0), (8) 

where are the parameters which obey the Bethe ansatz equations. We consider the wave 
function (8) in the sector Xi < X2 < . . . < xm- 

Using the definition of the operators B[qi) in terms of the S- matrices, reordering the S- 
matrices in the product Hi B{qi), we rewrite the equation (8) in terms of the new monodromy 
matrices Tj = Sio-^^Sio^ . . . Sio^^, i = 1, . . . L, depending on the site i and acting in the new 
quantum space (Oi, . . . Om), as follows: 

ip{xi, . . . Xm) = {0\{xi . . . xm\TiT2 . . . Tl\0)l\{11...1}m), 
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where the second average corresponds to the space (Oi, . . . Om)- We evaluate the average in 
the quantum space 1, 2, ... L, use the symmetry 0-H-l,A-H-D,B-ir^C and transform the 
operators acting in the space (Oi, . . . Om) to the operators in the F-basis. One should also take 
into account that the action of the operator F to the vacuum is trivial: F|0) = |0). 

Thus we have to calculate the following average over the states on the lattice with M sites 
Oi,...Om: 

^{x,,...xm) = ({ll...l}M|A^(6)^^(6)-i^''(exJ...i^^(e.M)-^''(a)|0), (9) 

where the operators B^{^xJ, i = 1, . . . M are located at the positions xi, . . . xm and the oper- 
ators and correspond to the new transfer matrix of the form Tj = Sio^ Sic,.^ . . . Smj^j . The 
factorizing operator and the operators in the F-basis are also defined in the space (Oi, . . . Om)- 
In contrast with the similar expression with the operator A and B in the usual basis, one can 
obtain the compact expression for the average (9) with the operators in the F-basis. 

Since the matrix A^ is diagonal, the next step to obtain the wave function is to use the 
equation (5) to rewrite the average (9) as a sum over the permutations. In fact one obtains: 

^{x,,...Xm)= E ({ll...l}M|^^(6M''(e2)-i^Pl(e.J-i^PM(e.M)-^''(a)|0). (10) 

Now using eq.(6) we commute all the operators Bpi{^xi) to the left. The action of the operators 
A^{^i), I 7^ Xi to the right- hand state (to the vacuum) is known, so we obtain the following 
compact expression (in the sector Xi < x^ < ■ ■ ■ < xm)'- 

n n c(6 - qj)m-i]M\Bp,{u)...BpM{i.M- (n) 

l = l^l^Xi j = l 

Note that one have to change the signs of all the spectral parameters in the relation (6) and 
the formula for the action of the operators ^^(6) to the vacuum due to the definition of the 
new monodromy matrices acting in the new quantum space (Oi, . . . Om)- The last step is 
to evaluate the average at the end of the last equation. Acting consequently by the operators 
Bpi to the right we obtain the expression: 

^M{^x,q\P) = ({ll...l}M|5pi(ea:J-5pM(e.M)|0) 

= n k^x. - qpi) n - 5P.) n v (12) 

1=1 i>j i>j ^{'ipt 'iPj) 

The sum over the permutations P of this expression gives the partition function of the six- 
vertex model with domain-wall boundary conditions J2p^Mi^x,(l\P) = ^m{C,x,(i)- It is in- 
teresting to obtain the determinant representation for this function ^m{^x,(i) [10]i[11] start- 
ing from the representation (12). In fact, representing the sum over the permutations in 
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J2p^Mi^x,<l\P) as J2i=il2p:PM=i^M{Cx,(l\P), and considering separately the dependence on 
the variables ^^'^ 1pm = Qi, we obtain the following recurrence relation for the function 

*M({e.}, {q}) = EKW - <l^) n '!5r^V ^M-l({e.Ja^M, {q^}m), 
i=l ajti '"vyj 

which coincides with the well known recurrence relation which determines the determinant 
expression for ^m{^x: q) (for example, see Appendix B of ref.[8]). 

Substituting the expression (12) into the equation (11) and performing the cancellations 
of similar terms we easily obtain the following result: 

^{Xi, . . . Xm) = A{P)(f)pi{xi)(j)p2{x2) . . . (pPMixM), (13) 
P 

where the functions (f)pi{xi) and the amplitude A{P) are equal to 

L Ml 

<f>Pi{xi) = n 5(6 - qPiMx, - qpi), A{P) = n ^( ^ • (14) 

The expression (13) for the wave function is the final result of the present paper. In order to 
compare this result with the results obtained long time ago by the other authors it is useful 
to represent the function (f)pi{xi) in the following form: 

M^i) = f n 5(6 - qpi)) {c-\^x, - qpMx, - qpi)) n 5"'(6 - qpi)- (i5) 

\l = l / l<Xi 

Thus, up to the normalization factor the wave function (13) coincides with the wave function 
obtained by the other authors (for example, sec [3], [1]). 

Let us obtain the Bcthc ansatz equations for the inhomogeneous six-vertex model starting 
from the wave function (13). We have to continue the wave function out of the interval {l,L) 
and impose the periodic boundary conditions of the form: 

ij{xi, . . . Xm) = i^ixi + L,X2,... Xm) = ^{x2-, ■ ■ ■ xm, Xi + L) (16) 

in the sector xi < X2 < ■ ■ ■ < xm- Substituting the wave function (13) into the equation (16) 
we obtain: 

YA{P)(t)pi{xi) . ..^pm{xm) = Xl^(^C')^Pi(a^i + L)(j)p2{x2) . ..(Ppm{xm), 
p p 

where C - is the cyclic permutation (CI = 2, C2 — 3. ...CM — 1). This equation gives the 
following condition for the amplitude A{P) (14): 

A{P)/A{PC)^l[c-\^,-qp^). (17) 
1=1 
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One can easily verify that the equations (17) are equivalent to the usual Bethe ansatz equations 
for the six- vertex model with the inhomogeneity parameters 

In conclusion, we derived the coordinate space wave function for the inhomogeneous six- 
vertex model from the Algebraic Bethe Ansatz. Our result is in agreement with the result 
first obtained by Yang and Gaudin for the eigenstate of the transfer matrix of the six-vertex 
model in the rational case. 
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